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A number of naturally occurring, proton-rich nuclides (the p-nuclei) cannot be made in the s- and 
r-process. It has been found that massive stars can produce p-nuclei through photodisintegration 
of pre-existing intermediate and heavy nuclei. This so-called y-process requires sufficiently high 
temperatures and occurs in pre-explosive or explosive O/Ne burning, depending on the mass of 
the star. Although the y-process has been successful in producing a large range of p-nuclei, two 
mass regions remain problematic, A < 1 10 and 150 < A < 165, where a number of p-nuclei are 
severely underproduced. The origin of the problems is yet to be identified. A large number 
of unstable nuclei with only theoretically predicted reaction rates are included in the reaction 
network and thus the nuclear input may involve uncertainties. Deficiencies in charged-particle 
optical potentials at y-process temperatures have been found for nuclei at stability. On the other 
hand, the y-process conditions (temperature profiles, entropy of the O shell, seed composition) 
also sensitively depend on details of the stellar structure and evolution, as well as on the initial 
metallicity. Nevertheless, especially the deficient low-mass p-nuclei may call for an additional 
production process or site, such as production in (subChandrasekhar) type la supernovae. Also 
the (im)possibility of a synthesis in the rp- and vp-processes is discussed. Were this the case, 
the production of p-nuclei would be realized as a superposition of several different processes, not 
necessarily in the same site. 
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1. The case of the missing nuclides 

When B 2 FH [1] and A. G. W. Cameron [2] postulated the s- and r-processes for the production of 
intermediate and heavy nuclei beyond iron, they realized that a number of proton-rich isotopes can 
never be synthesized through sequences of only neutron captures and j3~ decays. This required 
the postulation of a third process. It was termed p-process because it was thought to proceed via 
proton captures at high temperature, perhaps even reaching (partial) (p,/)-(y,p) equilibrium. B 2 FH 
tentatively placed this nucleosynthesis process in the H-rich envelope of type II supernovae but it 
was later realized that the required temperatures are not acquired there [3,4]. This also shed doubts 
on the feasibility to use proton captures for producing all of the nuclides missing from the s- and 
r-process production. 

It is somewhat confusing that in the literature the name "p-process" is sometimes used for 
a proton capture process in the spirit of B 2 FH but also sometimes taken as a token subsuming 
whatever production mechanism(s) is/are found to be responsible for the p-nuclides. For an easier 
distinction of the production processes, perhaps modern nomenclature better focuses on naming 
the nuclides in question the p-nuclides (they were called "excluded isotopes" by [2]) and use ap- 
propriate names for the processes involved. 

Compared to the bulk of natural isotopes, the p-nuclei generally show abundances which are 
1 — 2 orders of magnitude lower. Historically there were 35 p-nuclides identified, with 74 Se being 
the lightest and 196 Hg the heaviest. It is to be noted, however, that this assignment depends on 
the state-of-the-art of the s-process models (just like the "observed" r-abundances depend on them) 
and also on estimates of r-process contributions (e.g., to 113 In and 115 Sn) [5,6]. Already B 2 FH 
realized that, e.g., the abundance of 164 Er stands out among its neighboring p-nuclides and may 
thus contain considerable contributions from the s-process. It was indeed found that there are large 
s-process contributions to 164 Er and also to 152 Gd [7], thus possibly removing them from the list 
of p-isotopes. If also the abundances of 113 In and 115 Sn can be explained by modifications of the 
s-process and/or contributions from the r-process [6], this would leave only 31 p-isotopes to be 
explained by other processes. 

A special obstacle not found in the investigation of the s- and r-process is the fact that there 
are no elements dominated by p-isotopes. Therefore, our knowledge of p-abundances are limited 
to solar system abundances derived from meteoritic material and terrestrial isotopic compositions. 
Before astronomical observations of isotopic abundances at the required discrimination level be- 
come feasible (if ever), it is impossible to determine p-abundances in stars of different metallicities 
and thus to obtain the galactic chemical evolution picture directly. On the other hand, depending on 
the actual p-production mechanism this may also be problematic for determining early s-process 
contributions. If the p-nuclides (or some of them) turn out to be primary (i.e., independent of metal- 
licity) or have a different dependence of production on metallicity than the s-process (perhaps by 
initially originating from r-nuclei), they may give a larger contribution to elemental abundances in 
old stars than the s-process. 

2. Processes and sites possibly contributing to the production of p-nuclei 

There are several possibilities to get to the proton-rich side. Sequences of proton captures may 
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reach a p-isotope from elements with lower charge number. They are suppressed by the Coulomb 
barriers, however, and it is not possible to arbitrarily compensate for that by just requiring higher 
plasma temperatures. At high temperature (y,p) reactions become faster than proton captures and 
prevent the build-up of proton-rich nuclides. Photodisintegrations are an alternative way to make 
p-nuclei, either by directly producing them through destruction of their neutron-richer isotopes 
through sequences of (y,n) reactions (these are the predominant photodisintegration processes for 
most stable nuclei), or by flows from heavier elements via (y,p) and (y,a) reactions and j3-decays. 
The latter possibility requires higher temperatures because unstable, proton-rich nuclei have to be 
produced through sequences of (y,n) reactions first. 

There are three ingredients influencing the resulting p-abundances and these three are differ- 
ently combined in the various sites proposed as the birthplace of the p-nuclides. The first one is, 
obviously, the temperature variation as a function of time, defining the timescale of the process 
and the peak temperature. This already points us to explosive conditions which accommodate both 
the necessary temperatures for photodisintegrations (or proton captures on highly charged nuclei) 
and a short timescale. The latter is required because it has to be avoided that too much material is 
transformed in order to achieve the tiny solar relative p-abundances (assuming that these are typi- 
cal). The second parameter is the proton density. While photodisintegrations and j3 decays are not 
sensitive to the proton abundance, proton captures are. With a high number of protons available, 
proton captures can prevail over (y,p) reactions even at high temperature. Last but nor least, we 
can vary the seed abundances, i.e. the number and composition of nuclei on which the photodis- 
integrations or proton captures act initially. If the final p-abundances depend sensitively on these 
seeds (and they do in most suggested sites), this implies that the p-nuclei are secondary and thus 
depending on some s- and/or r-process nuclides being already present in the material. 

Here, only the most relevant suggested sites and processes for the production of p-nuclides can 
be mentioned (for further details and additional alternatives, see [4] and references therein). The 
most favored process so far is the y-process occurring during explosive O/Ne-shell burning in mas- 
sive stars [8-1 1]. It was realized early that the abundances of most p-nuclei are inversely correlated 
with their photodisintegration rates [2, 8], pointing to an important contribution of photodisinte- 
gration. Massive stars provide the required conditions of transforming s- and r-process material al- 
ready present in the protostellar cloud. During the final core-collapse supernova a Shockwave ejects 
and heats the outer layers of the star, in just the right amount needed to produce p-nuclei through 
photodisintegration. (Stars with higher mass may reach y-process temperatures of 2 < T$ < 3 al- 
ready pre-explosively [11] although some of those pre-explosive p-nuclei may be destroyed again 
in the explosion.) The y-process occurs naturally in massive stars and does not require any artifical 
fine-tuning. It is appealing because it can co-produce the bulk of p-nuclides within a single site. 
All p-nuclei are secondary. The very rare 138 La can be produced through neutrino reactions with 
neutrinos emitted by the nascent neutron star emerging from the core collapse (the v-process) [12]. 
The equally rare 180m Ta probably also received a large contribution from the v-process. Neverthe- 
less, there are longstanding shortcomings in the production of certain mass ranges which persist 
also in current models. The light p-nuclei with A < 110 are strongly underproduced and also the 
mass range 150 < A < 165 seems problematic. This has triggered a number of investigations in 
astrophysics and nuclear physics aimed at resolving these deficiences. 

An interesting alternative is the possible production of p-nuclei in hot, proton-rich environ- 



3 



Origin of the p-Nuclei 



Thomas Rauscher 



ments. At least the light p-nuclides could then be made by proton captures. In combination with 
a modified seed abundance (strong enrichment in s-process nuclei), again the whole range of p- 
nuclei could be synthesized possibly within one appropriate site. Such an environment is provided 
in the thermonuclear explosion of mass-accreting white dwarfs [13, 14]. Exploratory, parameter- 
ized calculations for the canonical type la supernova model - the explosion of a white dwarf (WD) 
after it has accreted enough mass from a companion star to reach the Chandrasekhar limit - also 
found an underproduction of light p-nuclei as well as of 138 La and 180m Ta, even when assuming 
a seed enrichment of factors 10 3 — 10 4 in s-process nuclei from an AGB companion [4]. The 
synthesis of p-nuclei proceeds via a /-process also in these types of models. Another alternative 
is a subclass of type la supernovae which is supposed to be caused by the disruption of a sub- 
Chandrasekhar WD (mainly composed of C and O) due to a thermonuclear runaway in a He-rich 
accretion layer [15, 16]. High neutron fluxes are built up in the early phase of the explosion and 
a weak r-process ensues. Once the temperature exceeds T$ ps 2 photodisintegrations take over and 
move the nucleosynthesis to the proton-rich side where two processes act: the y-process, as de- 
scribed above, and additionally rapid proton captures on unstable nuclei at 3 < Tg < 3.5. The latter 
is somewhat similar to the rp-process but at much lower proton densities and thus closer to stabil- 
ity. The proton captures are hindered by (y,p)-reactions but nuclei with low capture <2-value can 
be bypassed by (n,p) reactions due to the available free neutrons [4, 15]. This so-called pn-process 
can efficiently produce the light p-nuclides from Se to Ru but it overproduces them in relation to 
the heavier ones. Again, a strong increase in the photodisintegration seed abundances would be 
required to produce all p-nuclei at solar relative abundances. (It has to be noted that here the nuclei 
produced in the pn-process would be primary whereas the others are secondary.) Both WD scenar- 
ios (canonical and sub-Chandrasekhar mass type la supernovae) suffer from the fact that they are 
difficult to simulate and self-consistent hydrodynamical models including accretion, pre-explosive 
burning, explosion, and explosive burning in turbulent layers are missing. A complete model would 
also allow to study the actual amount of seed enhancement (if any) in these environments. Recent 
calculations based on trajectories from 2D simulations of WD explosions confirm the feasibility of 
p-nucleus production in these sites [16]. 

Considering even proton-richer conditions, two additional types of rapid proton capture pro- 
cesses may also yield p-nuclides from the decay of very proton-rich progenitor nuclei. Both of them 
imply that p-production is primary. The rp-process occurs in explosive H- and He-burning on the 
surface of a mass-accreting neutron star and involves sequences of proton captures and /3 -decays 
along the proton dripline [17]. How far the burning can move up the nuclear chart depends on 
details of the hydrodynamics (convection) and the amount of accreted protons. There is a definitive 
endpoint, however, when the rp-process path runs into the region of Te a-emitters [18]. Therefore, 
if the rp-process actually runs that far, only p-nuclides with A < 1 10 can be reached through decays 
of very proton-rich progenitors. Currently it is also unclear whether the produced nuclides can be 
ejected into the interstellar medium or whether they are trapped in the gravitational field and even- 
tual only modify the surface composition of the neutron star. On the other hand, very proton-rich 
conditions are found in the innermost ejected layers of a core-collapse supernova due to the inter- 
action of v e with the ejected, hot, and dense matter at early times. The hot matter freezes out from 
nuclear statistical equilibrium and sequences of proton captures and /3 -decays ensue, similar as in 
the late phase of the rp-process. The nucleosynthesis timescale is given by the explosion and subse- 
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quent freeze-out and is shorter than in the rp-process. The flow towards heavier elements would be 
hindered by nuclides with low proton capture Q- value and long /3 -decay halflives, similar as in the 
pn-process. Again, (n,p) reactions accelerate the upward flow. The neutrons stem from the reaction 
v e + P — > n + e + . Although at early times the v e flux is small, the sheer number of free protons 
guarantees a constant neutron flux. Thus, this process was termed the vp-process [19,20]. Like 
the rp- and pn-processes it may contribute to the light p-nuclides via decays of very proton-rich 
species. The details depend sensitively on the explosion mechanism, the neutrino emission from 
the proto-neutron star, and the hydrodynamics governing the ejecta motion. 

Historically, the longlived radionuclides 92 Nb, 97 Tc, 98 Tc, 146 Sm are not labelled as p-nuclides 
but they cannot be produced in the s- and r-process, either. Measurements in meteoritic material 
firmly established the presence of 92 Nb and 146 Sm at the formation of the solar system whereas 
only upper limits can be specified for the Tc isotopes [21]. The measured abundances of the 
above radionuclides are consistent within galactic chemical evolution models (i.e., they yield a 
unique chemical evolution history). The fact that 92 Nb existed at the formation of the solar system 
[22, 23] puts tight constraints on the contributions of the rp-, np-, and vp-processes to the light 
p-nuclei. The p-radionuclide 92 Nb cannot be reached by decays from the proton-rich side because 
it is shielded by its stable isobar 92 Mo. Combining the abundance information on the stable p- 
nuclei with those of the p-radionuclides the 92 Nb/ 92 Mo production ratio rNbMo in the Galaxy can 
be inferred: 1.7 ±0.5 x 10~ 3 < rNbMo < 0.1 [21]. The 92 Mo production in the /-process is too low 
by a factor of about 10. If the remaining 90% stem from decays of very proton-rich nuclei (which 
cannot produce 92 Nb) this would yield a production ratio of 1.5 ± 1.0 x 10~ 4 which is clearly 
outside the permitted range for rNbMo- This puts severe doubts on any significant contribution 
of the rapid proton capture processes to nuclei with A > 92. (A similar argument applies to the 
production of 98 Tc but it is at higher mass and its early solar system abundance is not well known.) 

3. Nuclear physics uncertainties 

Despite of the number of suggested sites, the nuclear physics underlying the p-production is similar 
in all of them, especially if a contribution of rapid proton capture processes far from stability can 
be ruled out. The nuclei involved in the /-process are the stable nuclides and moderately unstable, 
proton-rich nuclei. The masses have been measured and thus the reaction Q- values are well known. 
Halflives are known in principle but the electron captures and j3 + -decays need to be modified in the 
stellar plasma by the application of theoretical corrections for ionization and thermal excitation. 

Stellar correction factors act strongly also in the photodisintegration reactions (and captures). 
At the relatively high temperatures encountered in p-nucleosynthesis, a considerable fraction of nu- 
clei is thermally excited. Additional reactions proceeding from the excited states of nuclei become 
possible. Thus, laboratory measurements always only include a fraction of the astrophysically rel- 
evant transitions. To minimize the correction effects, reactions should always be determined in the 
direction of positive (?- value [24]. It can be shown, however, that also charged particle captures 
with negative Q-values show smaller corrections than their photodisintegration counterparts, due 
to the Coulomb suppression of the stellar enhancement [25]. Similar considerations apply for (p,n) 
and (n,p) reactions [26]. Unfortunately, measurements of the (p,y) and (a, 7) cross sections close to 
the astrophysically relevant energy windows of several MeV [27] at 2 < Tg < 3 (or 3 < Tg < 3.5 for 
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proton captures in the pn- and vp-processes) are scarce even for stable nuclei [28,29]. Therefore 
the full y-process reaction network includes a majority of rates based on theoretical predictions. 
Required are (n,y), (p,y), and (a,y) reactions and their inverses on stable nuclei and up to about 15 
units into the region of proton-rich, unstable nuclei (proton-richer nuclei are required in the rapid 
proton capture processes [30]). Not only the pn- and vp-processes contain (n,p) reactions, they are 
also important in the y-process where they are initiated by the neutrons released through (y,n) and 
speed up the matter flow back to stability [31]. It should be kept in mind that the theoretical correc- 
tions for the thermally excited nuclei have to be predicted even if a full database of experimental 
(ground state) rates were available. 

Charged particle captures and photodisintegrations are only sensitive to the optical potential at 
the astrophysically relevant low energies. These optical potentials are usually derived from elastic 
scattering at higher energies and are not well constrained around the Coulomb barrier. Especially 
the imaginary part should be energy-dependent (see [24] for more details on optical potentials 
in astrophysical rate predictions and the influence of other nuclear properties). Comparison of 
theoretical predictions with the few available data at low energy revealed a systematic 2-3 times 
overprediction of the (a,y) cross sections when using the "standard" potential of [35] (see, e.g., 
[28, 29, 36-39] and references therein). So far, the only known example of a larger deviation 
was found in 144 Sm(a,y) 148 Gd (determining the 144 Sm/ 146 Sm production ratio in the y-process) 
where the measured cross section is lower than the standard prediction by more than an order of 
magnitude at astrophysical energies [40]. Since this conclusion hinges only on a few datapoints 
with comparatively large errors, an independent (re-)measurement of this reaction is important. 

While the standard a+nucleus potential is a purely phenomenological one, the optical poten- 
tials for protons and neutrons used in the prediction of astrophysical rates and also the interpretation 
of nuclear data are based on a more microscopic treatment, a Bruckner-Hartree-Fock calculation 
with the Local Density Approximation including nuclear matter density distributions which are in 
turn derived from microscopic calculations [24,32-34]. Nevertheless, they also may show some 
deficiencies at low energies because also such descriptions use scattering data to constrain their 
free parameters [41]. In comparison to measured low energy (p,y) cross sections, these potentials 
fare much better than the one for the a+nucleus optical potential, being often in very good agree- 
ment with the data, sometimes being off by a maximal factor of two (see, e.g., [24,42,43] and 
references therein). In general, the uncertainty in the astrophysical rates caused by the nucleon 
optical potentials seems to be much lower than for the a-capture rates. Recent (p,y) and (p,n) 
data of higher precision close to the astrophysically relevant energy window, however, revealed a 
possible need for modification of the imaginary part of the nucleon potentials [24,43,44]. A con- 
sistently increased strength of the imaginary part at low energies improves the reproduction of the 
experimental data for a number of reactions. 

4. Deflections and branchings 

What is the actual impact of uncertainties in the photodisintegration rates on the production of the 
p-nuclides? To answer this question - which is also relevant for determining the main targets of 
future investigations - it has to be realized that not all reactions are equally important in all sections 
of the nucleosynthesis network. Furthermore, their relative importance changes with the plasma 
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Figure 1: Deflection points in the photodisintegration path for three stellar temperatures; natural isotopes 
are marked by squares. At low temperature, often the (y,n) reactions cease within an isotopic chain before 
reaching a path deflection. 

temperature reached in the respective mass zone where p-nuclei are produced. In a /-process, light 
p-nuclei are predominantly produced at higher temperatures (allowing efficient photodisintegra- 
tion of the nuclei around mass A « 100) whereas the production maximum of the heavy species 
lies towards the lower end of the temperature range [4,31]. Neutron captures and especially (y,n) 
reactions are important throughout the y-process network as the photodisintegration of stable nu- 
clides commences with (y,n) reactions until sufficiently proton-rich nuclei have been produced and 
(y,p) or (y,a) reactions become faster. The released neutrons can be captured again by other nuclei 
and push the reaction path back to stability in the region of the light p-nuclides. 

Studies of the y-process with systematic variations of the (y,p) and (y,a) rates by constant 
factors [31] or using a selection of different optical potentials for the (y,a) rates [45] confirm 
that uncertainties in the (y,p) rates mainly affect the lower half, whereas those in the (y,a) rates 
affect the upper half of the p-nuclide mass range. This can be easily understood by examination of 
the reaction paths and the competition between the photodisintegration reactions. In each isotopic 
chain we commence with initial (y,n) on stable isotopes and move towards the proton-rich side. The 
proton-richer a nucleus, the slower the (y,n) rate while (y,p) and/or (y,Of) rates increase. At a certain 
isotope within the isotopic chain, a charged particle emission rate will become faster than (y,n) and 
thus deflect the reaction sequence to lower charge number. Historically, this endpoint of the (y,n) 
chain has been called "branching", inspired by the branchings in the s-process path [8,45]. A more 
appropriate term would be "deflection (point)", though, because - unlike in the s-process - the 
reaction path does not split into two branches. The relative changes of the photodisintegration rates 
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Figure 2: Changes in the deflections when dividing the (y,0!) rate by a factor of 10. Only nuclei are marked 
which show different deflections than the standard case shown in Fig. 1. 



from one isotope to the next are so large that at each isotope only one of the emissions dominates. 
There are very few cases of two (or even three) types of emissions being comparable and those 
strongly depend on the optical potentials used. Examination of the deflection points easily shows 
that at higher mass we encounter (y,a) deflections whereas at lower mass most deflections are 
caused by (y,p) due to the distribution of reaction Q- values and Coulomb barriers [45]. 

Figure 1 shows the /-process deflection points for three relevant temperatures. The calcula- 
tions have been performed with the SMARAGD code [24,48] using the standard optical potentials 
of [32] and [35]. They supersede the results obtained in [45]. It is obvious that the higher the plasma 
temperature the further into the proton-rich side (y,n) reactions can act. At low temperature, even 
the (y,n) reaction sequence may not be able to move much beyond the stable isotopes because it 
becomes too slow compared to the explosive timescale (and neutron capture may be faster). Ob- 
viously, (y,p) and (y,a) are even slower in those cases and the photodisintegrations are "stuck", 
just slightly reordering stable abundances. For the plots presented here, photodisintegration rates 
slower than 0.2 s _1 were neglected. 

The above said, it has become obvious that not all possible reactions in the y-process network 
have to be known with high accuracy. Rather, only the dominant reaction sequence has to be known 
accurately, and within this sequence the slowest reaction as it determines the amount of processing 
within the given short timescale of the explosive process. A large error in a rate which is slower 
by several orders of magnitude than the dominating rate is inconsequential. Charged particle rates 
are important only at or close to the deflection points because they are many orders of magnitude 
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Figure 3: Uncertain deflections and branchings when simultaneously varying (y,n), (y,p), and (y,a) rates. 

smaller than the (y,n) rates on neutron-richer isotopes. The impact of suppressing all (y,ot) rates 
by a factor of ten is shown in Fig. 2. Only a comparatively small number of deflections are altered. 
Despite of the result of [40], such a large change may be exaggerated. Other reaction data suggested 
suppressions by factors of 2-3 (see previous section). When recalculating the deflections using the 
potential of [46,47] (which provides a good reproduction of the low-energy reaction data), there are 
even fewer modifications. Compared to what is shown in Fig. 1 an additional (y,n)/((y,p)+(y,a)) 
branching at 187 Ir appears at Tg = 2.5, and a (y,p)/((y,n)+(y,a)) branching at 100 Pd is replaced by a 
(y,p) deflection at = 3.0, thereby eliminating the (y,p) deflection at 98 Pd. 

A simultaneous variation of all photodisintegration rates within reasonably chosen limits, en- 
compassing the maximally expected rate uncertainty, can be taken as the extreme case. Figure 
3 shows the changed deflections (and branchings) when simultaneously varying the (y,n) rate by 
factors 1, 2, and 0.5, the (y,p) rate by factors 1, 2, and 1 /3, and the (y,a) rate by factors 1 and 0.1 
(for each nucleus separately). Again, the (y,n) photodisintegrations cannot efficiently destroy the 
stable isotopes of some elements at low temperature. 

5. Conclusion 

We have come a long way since [1,2] but the mystery of the origin of the p-nuclides is still with 
us. Although there are considerable uncertainties in the astrophysical (hydrodynamical) modeling 
of the sites possibly producing p-nuclei, a sound base of nuclear reaction rates is essential for all 
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such investigations. As long as an experimental determination of the rates around the deflection 
points is impossible, measurements of low energy cross sections of stable nuclides are essential to 
test and improve the theoretical calculations. This has been underlined by the recent results regard- 
ing optical potentials for the interaction of charged nuclei. Further measurements (at even lower 
energy) are highly desireable and, along with improved, self-consistent hydrodynamic simulations 
of the possible production sites, will gradually improve our understanding of p-nucleosynthesis. 
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